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We propose using a nonlinear phase shift interferometric converter (NPSIC) (a new 
Ph ! device) for lumped compensation of the nonlinearity in optical fibers. The NPSIC is 

a nonlinear analog of the Mach-Zehnder interferometer and provides a way to control 
the sign of the nonlinear phase shift. We investigate a potential use of NPSIC for 
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Abstract 



' compensation of the nonlinearity to make a dispersion-managed system closer to an 



ideal linear system. More importantly NPSIC can be used to essentially improve 
single channel capacity in the nonlinear regime. 



P : OCIS codes: 060.2330, 060.5530, 060.4370, 190.5530, 260.2030. 



The recent inventionlll and testingihi of the dispersion management technique demon- 
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^ I strated the effectiveness of this approach for high speed communications. Optical pulse dy- 
namics in fiber links with dispersion management are governed by the nonlinear Schrodinger 
equation (NLS) with periodic coefficients: 

iUz + d{z)utt + cr{z)\u\'^u = 0, (1) 

where z is the propagation distance, u is an optical pulse amplitude, d{z) = — 1/32(-2), 
^2(2) is a first order group-velocity dispersion, a{z) = {27m2)/{XoAeff{z)) is the nonlinear 
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coefficient, n2 is the nonlinear refractive index, Aq = 1-55 fim is the carrier wavelength, Ae/f 
is the effective fiber area that in general case depends on z. 

On short scales, the dispersion managed (DM) system is practically linear. Linear trans- 
mission in an optical fiber is limited by a nonlinear distance z < Zni = (cr|MoP)~^, which 
is determined by the Kerr nonlinearity a and the characteristic pulse power Iwop. The 
characteristic power cannot be choosen too small in order to maintain appropriate value of 
the signal-to-noise ratio. It is natural to attempt to extend the scale of the applicability 
of the linear regime. This can be achieved through use of a new optical fiber with lower 
Kerr nonlinearityi'0. Another obvious approach is "nonlinearity management", which was 
considered ini. However, the semiconductor material waveguides with negative Kerr non- 
linearity proposed as an element for compensation of nonlinear phase shift are currently not 
practical. In this letter we consider lumped compensation of the nonlinearity, the analog of 
lumped compensation of the chromatic dispersion by means of chirped fiber gratings. We 
suggest the use of a nonlinear phase shift interferometric converter (NPSIC) as presented in 
Fig.l. The NPSIC consists of a silica-based fiber 1, a highly nonlinear fiber 2 (which could 
be chalcogenide glass based, having a Kerr coefficient about 400 times or even much more 
than that for silica, see e.g. referencei), a linear amplifier A with amplitude amplification 
coefficient G, and two directional couplers 1,2. It is assumed that an optical terminator is 
installed at the end of the fiber 2 after coupler 2 to prevent beam reflection from the end of 
fiber 2. The incident light with amplitude Guq is split by the directional coupler 1 into two 
beams with amplitudes aiGuo and a2Guoe'^'^^'^ (ai > 0, a2 > 0), corresponding to the fibers 
1 and 2, respectively. We assume that total power is conserved: af + al = 1 but a subsequent 
consideration can be easily generalized to include the directional coupler's and fiber's losses. 
The extra phase 7r/2 in the second fiber is due to light spliting in the directional coupler 
(see e.g. reference^). The optical lengths /i, I2 in fibers 1 and 2 between couplers 1 and 2 
are chosen in such a way that they provide a zero phase difference at the input of coupler 2: 
^eff,ih — '"■e//,2^2 = 0, whcrc ^^e//,l , '"'e//,2 are effective linear refraction indexes in fibers 1,2. 
We assume that li, I2 are small enough so that we can neglect the influence of dispersion in 
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both fibers and the nonhnear phase shift in fiber 1. Thus the amplitudes of optical beams at 
the input of the coupler 2 are given by aiGuQ and a2Guoe'^'^^'^~^'^'^"^ , where the nonlinear phase 
shift (pni = <yhniO'2Cr'^\'^o\'^h and ahni correspond to the value of a in the highly nonlinear fiber 
2. Assuming that the propagation constants are the same for symmetric and antisymmetric 
modes of the coupler, the coupled-wave equation describing mode evolution in directional 
couplers can be written aJl^ 

{ui)z = inu2, {u2)z = inui, (2) 
where 

Ulout = Ul{Zout), (3) 

zq, Zout are the coordinates of coupler 2 input and output respectively; Uiout is the optical 
beam amplitude in fiber 1 at the exit of the coupler 2 and k is a coupling coefficient which is 
assumed to be real in the lossless model. The solution of (0), (§) shows that NPSIC converts 
the Mo into the output signal Uiout as follows 

uiout = Guo{ai cos 0c - a2e"^"' sin^J, (4) 

where 0c = i^izout — zq). We assume that 0„/ = crhni(i2G'^\uo\'^h ^ 1 and expand e***^"' in the 
Eq. (Q). Neglecting 0(0^^) terms results in 

Uiout = uiin[i- 7 . , ), (5) 

Oi COS 0c — a2 sm 0c 

where uun = Guo{ai cos (pc — 02 sin 0c) is the output amplitude Uiout of the linear system. 
Thus by changing the NPSIC parameters ai, 0c, /, G, one can control the sign and magnitude 
of the nonlinear phase shift. NPSIC can be considered as a nonlinear version of the Mach- 
Zehnder interferometer. As a typical example we take G = a/50, ai = 1/3, 0c = 0.2 which 
results in Uiout — Uo{1.0 — i60. cThniluol'^h)- Thus to completely compensate a nonlinear phase 
shift a\uo\'^Li in Li = 40 km line of a silica fiber we need to use NPSIC with I2 = 1.6m 



of a highly nonhnear fiber provided that ahni/o' = 400 (for this value of atni the nonlinear 
absorption is negligible in currently available chalcogenide glassed). 

This estimate is true if the power |up is constant throughout the propagation. In practice, 
power is not constant, due to fiber chromatic dispersion and losses. Therefore, nonlinearity 
compensators must be distributed along the fiber with a separation distance less than the 
dispersion length z^isp = r^/d and loss length zioss, where r is a typical pulse width. It is 
shown below that the effective length of the nonlinearity is increased as Zeff,ni ~ N'^Zni, 
where N is the number of lumped NPSIC on the dispersion map period. However, we 
demonstrate that for short pulses corresponding to strong dispersion management this ap- 
proach is efficient only at relatively large value of A^. Nevertheless, this method displays 
good performance in the nonlinear regime, as, for example, inserting only two compensating 
elements on the period of the dispersion map could increase the bit rate per channel by a 
factor of two. 

Consider DM system with step- wise periodical dispersion variation: d{z) = + d{z), 
d{z) = di for < z + mL < Li (standard monomode fiber) and d{z) = ^2 for Li < z + mL < 
Li + L2 (dispersion compensating fiber). Here do is the path-averaged dispersion, di, ^2 are 
the amplitudes of dispersion variation subject to a condition diLi + (i2-Z^2 = 0, L = Li + L2 
is a dispersion compensation period and m is an arbitrary integer number. a{z) = ai for 
< z + nL < Li and a{z) = 02 for Li < z + riL < Li + L2 corresponding to standard 
monomode and dispersion compensating fiber. We suppose that NPSIC units are located 
inside and at the ends of a dispersion compensated fiber at points z„ = mL + Li + nL2/N, 
n = 0, . . . , N, zq = mL + Li < Zi < . . . < z^ = {m + 1)L. At these points the value of 
u{t, z) experiences jump according to the Eq. (||): 

^(^' ^) L=.^.n+0 = - ^^e//,n \U\^U) L=,^^^_o, (6) 

where Zm,n = rnL + Li+ nL2/N; Zm,n — and Zm,n + mean the coordinate value just before 
and after the jump, respectively; the parameters ai, 0^, ^ G are chosen in such a way in 
order to to provide uun = u{t, z)\^=^^^^^o and (ThnialG^k sin (pc = ^e//,n = {criLi + a2L2)/N 
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for n=l,. . . , N-1 and -Ze/z.o = Zeff,N = + o"2-^2)/(2A^). Here the term — (criLi + (T2-L2) 

provides the compensation of the nonhnear phase shift both in standard monomode and 
dispersion compensating fibers. 

Assuming that the nonhnearity is small Zni ^ Zdisp, Zni ^ L, where L is a dispersion 
map period, one can express u in Fourier domain as a product of an exact solution of linear 
part of the Eq. (|l]), corresponding to mean- free dispersion d{z), on a slow function i^{u!,z): 
u{u,z) = i(j{uj,z)e Lq'^'-^'^'^^ (^ggg referenceS), where u{uj,z) = J^^u(t, z)e^'^''dt. is a 
slow function of z on a scale L which allows to integrate the Eq. (|l]) over the period L 
neglecting the slow dependence of on z and we get by taking into account the jumps in 



itjj{uj, {m + 1)L) — iipiuj, mL) — Llu d^ipif^^ rriL) 

+ 777^ / 4^{^^i,mL)il!{uj2,mL)ilj*{u3,mL) (7) 



(27r)2 . 

xKtot{^)S{^i + UJ2 — — uj^)dujiduj2duj2, = 0, 

where s = diL^, A = uj + - - cul and Ktot(A) = Ki{A) + K2{A). A kernel Ki{A) 
is equal to the usual kernel of the path-averaged equation0: 

Ki{A) = ai sinx/x, x = sA/2, ai = aiLi + cr2L2, (8) 
while i^2(A) accounts for the jumps in (^: 

K2{A) = -aL [sin x cot {x/N)]/N. (9) 

Suppose that Full Width at Half Maximum (FWHM) r of 1^(0, t) is subject to the condition 
s/{2Nt'^) <^ 1 then Ktot{A) can be rewritten as: 

Ktot{A) = ai(sinx/x)[xV(3Ar2) + 0(xVAr2)]. (10) 



Thus comparing (|TOD with the kernel (H) of the usual path-averaged equation we get an 



extra small factor (f^) /3 ~ L"^ / {NzdispY provided that is big enough to ensure the 



condition j^^^ <^ 1. For short pulses A^ should be a big number to make the system 



close to the linear one. Figure 2 shows the dependence of the FWHM Tout obtained af- 
ter propagation of the initial zero-chirp Gaussian pulse with rj„j = 10 over a typical 
transoceanic distance km on the number of NPSIC units + 1. This dependence is 
obtained by numerical intergation of the NLS equation (|T]). The pulse is lunched to the 
DM fiber at z = Li/2 with the peak power |mP = 2mW. DM system parameters are 
do = 0, di = lO.Ops'^/km, d2 = — 101.6 ps"^/ km, Li = 40 km, L2 = —diLi/d2, ci = 
0.0013 {km mW)~^ , (T2 = 0.00405 (/cm mW^)~^. It is seen that the system can be approx- 
imately considered linear for ~ 30 in which case the nonlinear correction to Tout is less 
5%. 

A strong variation of Tout in Fig- 2 for small A^ is due to the attraction of the initial 
Gaussian pulse to soliton solution ipi^y ™L) = ipQ{u!)e'^"^^^ of path-averaged Eq. (0). Here A 
is a soliton propagation constant. We refer to this soliton solution as the modified DM soliton 
(MDM soliton) by analogy with the DM soliton, that corresponds to the usual path-average 
equation (without nonlinearity management). For K2{A) = we recover the usual DM 
soliton of the path-averaged equationS where the DM soliton width tdm depends on s only 
for do (see e.g. reference^). For the above mentioned system parameters tdm — 21ps. 
Any shorter pulses experience strong distortion because of the nonlinearity for K2{A) = 0. 
Curve 2 in Fig. 3 shows the pulse power distribution after lO^fcm of propagation of the 
initial Gaussian pulse with rj„j = lOps (curve 1). Using NPSIC units we can control K2{A) 
and thus change tmdm of MDM soliton. Curve 3 shows the pulse power distribution after 
10^ km for a system with two NPSIC for each DM period L located at the ends of the 
dispersion compensating fiber with -Ze/z.o = -2e//,i = /(o^i-^i + o"2-^2)/2, where f=3/2. One 
can see that curve 3 is close both to the initial Gaussian pulse (curve 1) and to the MDM 
soliton (curve 4). The MDM soliton was obtained by numerical iteration of the equation (|^ 
with %lj{ijj,mL) = iljQ{uj)e'^'^^^, X = 0.00028 fcm~^. The numerical iteration scheme is similar 
to the one used in reference0 for usual path-averaged equation. By changing the factor 
/ one can control t^dm of the MDM soliton. Note that for / > 1, which corresponds to 
a negative average nonlinearity of the total system, the numerical iteration scheme finally 
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diverges indicating that the MDM sohton is a rather long-hved quasi-stable structure which 
however serves as an attractor of the pulse dynamics on a scale ~ 10^ km. 

In summary, the use of NPSIC elements would allow one to construct a nearly linear 
fiber transmission system. However this system requires many NPSIC elements. On the 
other hand the use of NPSIC elements offers the better advantages to fiber hnks operating 
in nonlinear regime. In particular, only a few elements per dipersion map period could 
dramatically reduce pulse width and potentially increase bit-rate. 

The authors thanks M. Chertkov and E.A. Kuznetsov for helpful discussions. 
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Figure captions: 



Fig.l. Schematic of the nonhnear phase shift interferometric converter. 

Fig. 2. FWHM Tout after pulse propagation over 10*^ km versus a number A'^ + 1 of NPSIC 

units. 

Fig. 3. Power distributions: an initial Gaussian pulse (curve 1); result of pulse propagation 
over 10^ km in a DM system with no NPSIC (curve 2); result of pulse propagation in a DM 
system with two NPSIC units for each period L (curve 3); the MDM sohton (curve 4). 



9 



FIGURES 



Fiber 2 




Figure 1 

I.R. Gabitovand RM. Lushnikov 
"Nonlinearity Management../' 

Subnnittecl to Optics Letters 



10 




N+1, Number of NPSIC 



Figure 2 

I.R. Gabitov and P. M. Lushnikov "Nonolinearity Management ..." 

Submitted to Opt. Lett. 



11 




time t, ps 



Figure 3. 

I.R. Gabitov and P.M. Lushnikov "Nonlinearity Management ..." 
Submitted to Optics Letters 



12 



